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Application  of  inverse  theory  to  the  observation  of  high-degree  five-minute  solar 
oscillations  has  led  to  the  detection  of  horizontal  flows  below  the  solar  surface  that 
are  a  combination  of  solar  rotation  and  giant  convection  cells.  The  distinctive  dis¬ 
placements  in  the  centroids  of  the  ridges  evident  in  the  power  diagrams  of  the  oscilla¬ 
tions  from  one  observing  day  to  the  next  arise  from  different  patterns  of  giant^cells 
being  rotated  into  view.  Such  observation  of  frequency  splittings  for  the  highfdegree 
oscillation  modes,  combined  with  refinements  in  the  inversion  of  the  data  using  optimal 
averaging  and  spectral  expansions,  has  shown  that  helioseismology  should  permit  detailed 
mapping  of  velocity  and  thermal  structures  below  the  solar  surface.  Extensive  theoreti¬ 
cal  studies  of  fully  compressible  magnetoconvection  have  shown  that  flows  are  indeed 
able  to  concentrate  magnetic  fields  into  concentrated  flux  sheets  that  are  substantially 
evacuated  of  gas.  The  magnetic  buoyancy  instabilities  have  been  extensively  studied.<^L_ 
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A.  OVERVIEW  OP  RESEARCH  GOALS 


The  research  supported  by  this  contract  has  been  directed 
-toward  the  study  of  the  nonlinear  dynamics  of  turbulent  convection 
below  the  surface  of  the  Sun,  and  the  ability  of  such  convection 
to  produce  strong  magnetic  fields  through  dynamo  action.  The 
coupling  of  convection  with  magnetic  fields  influences  much  of 
what  is  observed  in  the  solar  atmosphere,  and  so  too  most  of  the 
variations  seen  in  the  solar  wind.  Probing  of  the  underlying 
convection  is  essential  for  understanding  solar  variability. 

The  dynamics  of  the  convection  zone  has  been  studied  under 
this  contract  by  a  combination  of  observations  and  theory.  The 
five-minute  oscillations  of  the  Sun  have  been  used  to  probe 
velocity  and  thermal  structures  below  the  solar  surface.  This  has 
been  possible  because  the  presence  of  giant  convection  cells,  as 
predicted  by  theory,  leads  to  changes  in  the  apparent  frequencies 
of  the  acoustic  waves  responsible  for  the  oscillations.  Such 
frequency  shifts  with  time  have  been  measured  from  full-disk 
Doppler  observations  of  the  Sun  carried  out  from  Sac  Peak.  Appli¬ 
cation  of  inverse  theory  to  such  observed  frequency  splittings  of 
the  five-minute  oscillations  has  allowed  us  to  begin  to  map  the 
velocity  fields  with  depth  below  the  solar  surface.  Such  work 
in  helioseismology  has  been  complemented  by  nonlinear  models  of 
solar  convection  which  take  into  account  the  compressibility  of 
the  medium  and  permit  the  flows  to  penetrate  into  the  stably- 
stratified  atmosphere.  The  theoretical  models  have  also  studied 
the  interaction  of  convection  with  magnetic  fields,  dealing 
directly  for  the  first  time  with  the  possible  evacuation  of  the 
plasma  within  the  magnetic  flux  sheets  that  can  be  formed  by  the 
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convection.  The  theoretical  studies  have  proved  to  be  very 
useful  in  also  interpreting  the  observations  with  high  spatial 
resolution  of  both  granular  and  mesogranular  flows  on  the  solar 
surface  as  observed  from  both  Sac  Peak  and  from  the  Spacelab  2 
SOUP  instrument. 

B.  RESEARCH  ACCOMPLISHMENTS 

The  research  work  supported  in  substantial  part  by  this 
contract  has  resulted  in  30  published  (or  in  press)  papers  in 
refereed  scientific  journals  or  proceedings.  The  cover  pages 
with  abstracts  and  full  citations  are  attached  as  Papers  A  to  DD. 

A  further  16  published  abstracts  of  talks  at  national  conferences 
have  also  been  cited  in  the  quarterly  reports,  and  there  have  been 
over  28  other  seminars  given  describing  aspects  of  the  research. 

Thus  the  contract  has  supported  research  that  has  been  highly 
visible  and  thoroughly  documented  in  the  open  scientific  litera¬ 
ture. 

Group  1.  (Papers  A  to  E) :  Major  reviews  and  summary  discus¬ 
sions  of  helioseismology.  The  subject  of  using  the  five-minute 
oscillations  to  probe  the  interior  structure  and  dynamics  of  the 
Sun  has  become  an  area  of  intense  research  interest  and  enthusiasm 
in  solar  physics.  Helioseismology  has  now  also  been  embraced  by 
a  wide  group  of  physicists  as  a  topic  that  should  be  strongly 
encouraged  because  of  its  substantial  promise  in  permitting  study 
of  the  inside  of  a  star.  This  has  permitted  the  Global  Oscillations 
Network  Group  (GONG)  to  be  initiated  as  a  national  project  centered 
at  the  National  Solar  Observatory,  involving  both  Tucson  and  Sac 


Peak.  With  partial  assistance  from  this  contract,  we  have  been 
able  to  consider  helioseismology  from  a  broad  perspective  and 
to  help  explain  its  subtleties  in  a  number  of  prominent  review 
articles  within  this  group  of  papers.  In  particular,  Papers  C 
and  D  were  color  cover  articles  in  Science  and  in  Scientific 
American,  and  both  elicited  considerable  interest  and  discussion 
from  scientists  and  from  people  involved  with  national  policy. 

Group  2.  (Papers  F  to  H);  Search  for  solar  giant  cells 
and  study  of  rotation  with  depth.  Application  of  inverse  theory 
to  our  observations  of  high-degree  five-minute  oscillations  has 
permitted  us  to  detect  horizontal  flows  below  the  solar  surface 
that  are  a  combination  of  solar  rotation  and  giant  convective 
cells.  The  distinctive  displacements  in  the  centroids  of  the 
ridges  evident  in  the  power  diagrams  of  the  oscillations  from  one 
observing  day  to  the  next  arise  from  different  patterns  of  giant 
cells  being  rotated  into  view.  Our  observation  of  such  shifts  in 
the  frequencies  of  the  oscillations,  combined  with  our  refinements 
of  inversion  of  the  data  by  optimal-averaging  and  spectral  expan¬ 
sions,  is  serving  to  develop  helioseismology  as  a  detailed  approach 
for  studying  the  dynamics  of  the  solar  convection  zone. 


Group  3.  (Papers  N  to  0) :  Studies  of  active  regions  and 
flares.  These  two  papers  have  dealt  with  regions  possessing 
strong  magnetic  fields.  The  first,  involving  coordinated  SMM 
satellite  and  ground-based  observations,  revealed  that  active 
regions  possess  a  broad  spectrum  of  rapid  brightening  events, 
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suggesting  that  magnetic  field  reconnection  may  be  proceeding 
almost  continuously  and  stochastically.  The  second  paper  is 
unusual  in  describing  helioseismologic  observations  that  were 
carried  out  during  several  days  that  were  interrupted  by  a 
major  solar  flare.  This  permitted  the  first  evaluation  of  the 
flare  site  as  a  possible  source  of  acoustic  mode  excitation. 

Group  4.  (Papers  P  to  V) :  Theory  of  compressible  non¬ 
linear  convection  and  magnetoconvection.  Through  the  development 
of  compressible  theoretical  models  for  convection,  coupled  with 
the  use  of  supercomputers,  we  have  studied  the  interaction  of 
penetrative  convection  with  magnetic  fields.  In  particular,  we 
have  examined  how  convection  can  penetrate  downward  into  the 
stable  stratification  below  the  convection  zone,  finding  that 
lateral  swaying  of  the  convective  plumes  can  lead  to  the  excita¬ 
tion  of  gravity  waves  there.  That  region  is  thought  to  be  able 
to  support  magnetic  dynamo  action,  and  our  studies  have  provided 
estimates  of  the  ability  of  convection  to  supply  mechanical  energy 
into  that  region.  We  have  also  considered  the  behavior  of 
compressible  magnetoconvection,  finding  that  the  magnetic  fields 
are  concentrated  into  sheets  by  the  convective  motions.  These 
sheets  can  be  substantially  evacuated  of  gas  to  achieve  a  balance 
between  the  gas  pressure  outside  and  the  strong  magnetic  pressure 
within  the  sheet.  Remarkably  enough,  the  sheets  can  appear  almost 
as  stagnant  regions  of  fluid  surrounded  by  sheaths  of  fast  down¬ 
flow,  not  unlike  the  behavior  near  magnetic  pores  on  the  Sun. 

These  are  the  first  dynamical  calculations  to  consistently  treat 
compressibility  and  field  evacuation  in  magnetoconvection. 


Group  5.  (Papers  W  to  BB) :  Magnetic  buoyancy  instability. 


These  six  papers  study  the  nature  of  both  linear  and  nonlinear 
instability,  and  its  consequences,  on  magnetic  field  configura¬ 
tions  through  effects  of  magnetic  buoyancy.  This  research  has 
revealed  that  a  wide  range  of  magnetic  field  stratifications  can 
be  unstable,  for  there  are  about  three  varieties  of  instabilities 
that  can  disrupt  the  fields.  This  has  significant  implications 
upon  whether  it  is  possible  to  store  strong  magnetic  fields  at 
the  base  of  the  convection  zone  without  them  being  disrupted  by 
magnetic  buoyancy  instabilities.  These  studies  have  been  supported 
by  extensive  nonlinear  numerical  simulations  of  the  manner  in 
which  emerging  magnetic  flux  loops  twist  and  fragment. 

Group  6.  (Papers  CC  to  DP) ;  Nonlinear  dynamical  systems. 

These  two  papers  report  on  the  very  rich  behavior  observed  in 
doubly-dif fusive  convection  that  can  lead  to  periodic,  quasi- 
periodic  and  chaotic  motion.  The  first  paper  has  been  of  consid¬ 
erable  importance  in  revealing  that  partial  differential  equations 
can  allow  the  remarkable  transitions  previously  observed  with  much 
simpler  mathematical  systems.  This  has  bearing  on  the  forms  of 
time  dependence  that  should  be  achieved  with  magnetoconvection  as 
one  example.  The  second  paper  has  found  that  double-diffusive 
systems  can  admit  travelling  waves  for  the  convection,  and  has 
studied  their  relative  stability  to  standing  waves  normally 
studied  in  convective  systems. 


In  Solar  Seismology  from  Space  (eds.  R.K.  Ulrich  et  al.:  NASA/ 
Jet  Propulsion  Laboratory  Publ.  84-84),  pp.  7-39  (1984) 

PAPER  A. 

OTE1TIEW  OF  SOLAR  SEISMOLOCT;  OSCILLATIOMS  AS 
FIOIES  or  INTERNAL  STBOCTUBE  AND  PTNAH1CS  If  TIE  SOM 


JOIX  TOOM1E 

Department  of  Aatrophyslcal ,  Planetary  and  Atmospheric 
Science* ,  and  Joint  Institute  for  Laboratory  Aatrophytlcs , 
University  of  Colorado,  Boulder,  CO  80309,  U.S.A. 


ABSTRACT 

The  physical  nature  of  solar  oscillations  Is  reviewed.  The 
nosienclature  of  the  subject  and  the  techniques  used  to  Interpret  the 
oscillations  are  discussed.  Many  of  the  acoustic  and  gravity  waves 
that  can  be  observed  In  the  atnosphere  of  the  Sun  are  actually  resonant 
or  standing  nodes  of  the  Interior;  precise  neasureaents  of  the 
frequencies  of  such  nodes  allow  deductions  of  the  Internal  structure 
and  dynaalcs  of  this  atar.  The  scientific  objectives  of  such  studies 
of  solar  aelsalc  disturbances,  or  of  solar  seismology,  will  be  outlined. 
The  reasons  for  why  It  would  be  very  beneficial  to  carry  out  further 
observations  of  solar  oscillations  both  froa  ground-based  networks  and 
froa  space  will  be  discussed. 


In  Solar  Seismology  from  Space  (eds.  R.K.  Ulrich  et  al.:  NASA/ 

Jet  Propulsion  laboratory  Publ.  84-84),  pp.  49-78  (1984) 

PAPER  B. 

SOLAR  INVERSE  THEORY 
Douglas  Gough 

Institute  of  Astronomy  and  Departaent  of  Applied  Mathematics 
and  Theoretical  Physics,  University  of  Cambridge,  and  Joint 
Institute  for  Laboratory  Astrophysics,  University  of  Colorado 
and  National  Bureau  of  Standards,  Boulder,  Colorado  80309 


Helloselsmological  Inversion,  as  with  the  Inversion  of  any  other 
data,  is  divided  Into  three  phases.  The  first  Is  the  solution  of  the 
so-called  forward  problem:  namely,  the  calculation  of  the  elgenf requen- 
cles  of  a  theoretical  equilibrium  state.  The  second  is  an  attempt  to 
understand  the  results,  either  empirically  by  determining  how  those  fre 
quencles  vary  as  chosen  parameters  defining  the  equilibrium  model  are 
varied,  or  analytically  from  asymptotic  expansions  in  limiting  cases  of 
high  order  or  degree.  A  familiarity  with  at  least  the  qualitative  de¬ 
pendence  of  the  eigenfrequencles  on  various  properties  of  the  solar 
model  is  necessary  not  only  for  personal  enlightenment  but  also  for 
arming  oneself  to  Interpret  the  rather  more  abstract  third  phase.  That 
phase  Is  to  pose  and  solve  an  Inverse  problem,  which  seeks  to  find  a 
plausible  equilibrium  model  of  the  Sun  whose  eigenfrequencles  are  con¬ 
sistent  with  observation. 

The  three  phases  are  briefly  discussed  In  this  review,  and  the 
third,  which  Is  not  yet  widely  used  In  helloselsmology,  is  Illustrated 
with  some  selected  Inversions  of  artificial  solar  data. 
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6  September  1985,  Volume  229.  Number  4717 


Science.  229.  923-931 
(6  Sept.  1985) 


with  color  cover. 


PAPER  C. 


Seismology  of  the  Sun 

Jergen  Christensen-Dalsgaard.  Douglas  Gough.  Juri  Toomre 


J  Christensen  DaUgaard  i  lecturer  ,r  as'^ono 
n\)  »t  the  Astronomisk  Intutut  Aa'fcut  Lm*er*net 
DK  9000  Aarhus  C.  Denmark  and  a  rrvnrth  dip 
cute  at  NORDITa  Copenhagen  Denm*»a  D 
Cough  it  a  lecturer  in  asffonofrv  and  arr  r£  i^a'he 
matKt  at  the  Inv-tutc  of  Astronomy  a*';  the  Defa" 
mem  of  Applied  Mathematic*  and  Theo'etica  Rrw  v 
ic*  tni^crait)  of  Cambridge.  Cambridge  CB3 
OHA.  England  J  Toomre  it  a  profcivor  ir,  the 
Department  of  Astrophy  steal  Planetary  and  Aime- 
spheric  Science*  and  a  fello»  of  the  Joint  Intvjie 
for  Laboratory  Astrophysics  Ln.*erm>  of  Coiora 
do.  Boulder  Colorado  hGWr 


Summary  Oscillations  ot  th#  sun  make  it  possible  to  probe  the  inside  o<  •  star  The 
frequencies  of  the  oscillations  have  already  provided  measures  of  the  sound  speed 
and  the  rate  of  rotation  throughout  much  of  tha  sola'  interior.  These  quantities  are 
important  for  understanding  the  dynamics  of  the  magnetic  cycle  and  have  a  bearing 
on  tasting  general  relativity  by  planetary  precession  The  oscillation  frequencies  yield 
a  helium  abundance  that  is  consistent  with  cosmology,  but  they  reinforce  the  seventy 
of  the  neutnno  problem  They  should  soon  provide  an  important  standard  by  which  to 
calibrate  the  theory  ot  stellar  evolution 


Scientific  American.  253.  No.  3,  48-57, 
with  color  cover  (Sept.  1985) 

Helioseismology  PAPER  D, 

Acoustic  waves  within  the  sun  are  visible  as  oscillations 
on  the  solar  surface.  Their  pattern  and  period  hold  clues  to 
structure,  composition  and  dynamics  in  the  sun  s  interior 


by  John  W.  Leibacher,  Robert  W.  Noyes,  Juri  Toomre  and  Roger  K.  Ulrich 


OSCILLATIONS  OF  THE  SI' VS  SURFACE  in  the  rauH  of  tound  «i>b  monating 
la  Hi  interior,  Sere  four  of  Ike  10  millioo  retonincee  tkat  occur  In  tkt  tun  nrc  modeled  by 
computer  Surface  refiont  tkal  are  approaching  Ike  ok*er>er  are  colored  blur,  region!  tbit 
are  receding  are  colored  red.  la  actual  obtenation  tuck  turface  displacements  are  c>idtni 
at  Doppler  shifts  in  Ikr  wavelengths  of  light  abeorked  by  Ike  meting  gain  and  at  variation* 
la  krigktnma  The  ipatial  patlrm  and  tkc  ptriod  of  a  turface  otcillation  enable  investigators 
to  deduce  tbe  tkreo-dimentfonal  itrucrure  of  tbe  resonance  and  to  infer  properties  of  ibe 
tolar  Interior  For  each  otcillation  tbt  drgrrr  </)  and  Ibt  azimuthal  order  (ml  arr  given 
Degree  daertbes  the  pattern  in  lermt  of  Its  total  numbtr  of  eodrt  (circlet  along  »bicb  tbe 
aurfare  of  tbe  ipkere  b  motionless.  teea  krrr  at  handt  of  gray  bet- ten  tbe  lorn  of  color); 
azimuthal  order  it  aa  iadicatioa  of  the  lumber  of  eodet  that  ietenret  tbe  tolar  equator. 
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PAPER  E. 


in  Seismology  of  the  Sun  and  the  Distant  Stars 

(ed.  D.O.  Gough;  Reidel,  Dordrecht),  pp.  1-22  (1986) 


PROPERTIES  OF  SOLAR  OSCILLATIONS 


Juri  Toomre 

Joint  Institute  for  Laboratory  Astrophysics 

Department  of  Astrophysical,  Planetary  and  Atmospheric  Sciences 
University  of  Colorado 
Boulder,  Colorado  80309  USA 


ABSTRACT.  Many  of  the  oscillations  that  can  be  observed  in  the  atmosphere 
of  the  Sun  are  resonant  acoustic  or  gravity  modes  of  the  interior.  Accurate  mea- 
surement  of  the  frequencies  of  these  p  and  g  modes  permits  deductions  about  the 
internal  structure  and  dynamics  of  this  star.  Some  of  the  methods  of  interpretation, 
involving  a  close  interplay  between  observation  and  theory,  can  be  carried  over  to 
the  study  of  more  distant  stars. 

Mem.  Soc.  Astron,  I  tali  ana,  55,  153-161  (1984)  PAPER  F, 


ATTEMPT  TO  MEASURE  THE  SOLA*  SUBSURFACE  VELOCITY 


njjtx  dll 

National  Solar  Observatory,  Sunspot,  W  88349,  U.S.A.,  and 
Joint  Institute  for  Laboratory  Astrophysics,  University  of  Colorado, 

DOOCLAS  GOOCH 

Institute  of  Astronomy,  University  of  Cambridge,  Cambridge,  CB3  OKA, 
England,  and  Joint  Institute  for  Lsboratory  Astrophysics, 
University  of  Colorado 


JUKI  TOOK** 

Department  of  Astrophysical,  Planetary  and  Atmospheric  Sciences,  and 
Joint  Institute  for  Laboratory  Astrophysics,  University  of  Colorado 

Boulder,  CO  80309,  U.S.A. 


ABSTRACT  —  Five-minute  oscillation  modes  are  advected  by  horitontal 
velocities  below  the  solar  surface,  and  thus  can  be  used  as  probes 
of  rotation  and  large-scale  convective  flows.  Results  of  Inverse 
theory  applied  to  observations  of  hlgh-degree  modes  carried  out  on  alx 
separate  days  reveal  variations  In  horitontal  velocities  with  depth 
from  day  to  day  that  may  be  the  result  of  giant  convection  cells, 
though  noise  In  the  data  makes  this  Interpretation  somewhat  tentative. 


in  Solar  Seismology  from  Space  (eds.  R.K.  Ulrich  et  al.:  NASA/ 
Jet  Propulsion  Laboratory  Publ.  64-84),  pp.  79-93  (1984) 
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rotational  inversion  from  global  solar  oscillations 


Jsrgen  Chriscensen-Dalagsard 

High  Altitude  Observatory,  Boulder,  Colorado,  and 
NORDITA,  BlegdamcveJ  17,  Kebenhavn  0 

Douglas  Cough 

Institute  of  Astronomy  and  Department  of  Applied  Mathematics 
and  Theoretical  Physics,  University  of  Cambridge,  and  Joint 
Institute  for  Laboratory  Astrophysics,  University  of  Colorado 
and  National  Bureau  of  Standards,  Boulder,  Colorado 


We  investigate  the  degree  to  which  various  sets  of  solar  oscilla¬ 
tions  can  resolve  the  solar  internal  rotation.  Genuine  observations 
were  elmulated  by  the  following  procedure:  First  an  artificial  angular 
velocity  was  Invented  by  one  of  us,  and  froa  it  the  rotational  splitting 
of  a  aet  of  normal  node  a  ms  calculated;  to  that  ms  added  some  random 
noise.  The  result  was  treated  as  artificial  data  by  the  other  author, 
acting  as  an  observer,  who  attenpted  to  recover  the  rotation  law  by 
using  the  Beckus-Cllbert  optimal  averaging  procedure.  The  observer  knew 
neither  the  original  rotation  law  nor  the  mount  of  noise  that  had  been 
added.  Finally  his  conclusion  was  compared  with  the  actual  artificial 
angular  velocity. 

in  Solar  Seismology  from  Space  (eds.  R.K.  Ulrich  et  al.:  NASA/ 
Jet  Propulsion  Laboratory  )>ubl.  84-84),  pp.  95-111  (1984) 
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SENSITIVITY  OF  INFERRED  SOBPHOTOSPHERIC  VELOC ITT  FIELD  TO 
MODE  SELECTION,  ANALYSIS  TECBKIQPE  AND  NOISE 
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Joint  Institute  for  Laboratory  Astrophysics,  University  of 
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DOUGLAS  CODCB 

Institute  of  Astronomy,  and  Department  of  Applied  Mathematics 
and  Theoretical  Physics,  University  of  Cambridge,  Cambridge, 

CB3  OHA,  England,  and  Joint  Institute  for  Laboratory 
Astrophysics,  University  of  Colorado 

JOEI  TOONRE 

Department  of  Astrophys leal ,  Planetary  and  Atmospheric 
Sciences,  and  Joint  Institute  for  Laboratory  Astrophysics, 
University  of  Colorado,  Boulder,  CO  80309,  U.S.A. 


ABSTRACT 

The  horizontal  velocity  immediately  below  the  photosphere  esn  be 
Inferred  from  observetlons  of  hlgh-degree  solar  oscillations  by  an 
optimal-averaging  inversion  technique.  We  investigate  the  sensitivity 
of  the  results  to  various  details  of  both  the  inversion  end  the  deter¬ 
mination  of  the  frequencies.  The  results  are  shown  to  be  quite  stable 
to  the  choice  of  most  parameters,  suggesting  that  this  procedure  pro¬ 
duces  reliable  estimates  of  the  subsurface  velocity. 
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in  Solar  Seismology  from  Space  (eds.  R.K.  Ulrich  et  al. :  NASA/ 
Jet  Propulsion  Laboratory  Publ.  84-84),  pp.  271-277  (1984) 


THE  EFFECTS  OF  A  HEA1LT  1001  PIHT  CYCLE  OH  OBS  EE  VAT  IONS  OF  SOLAS 
OSCILLATIONS 

PAPER  I, 

Frank  Hill  A 

National  Solar  Observatory 
Sunspot,  N.M.,  and 

Joint  Institute  for  Laboratory  Astrophysics 
University  of  Colorado. 

Aicmff;  Power  spectra  of  window  functions  with  duty  cycles  between 
60S  and  99X  and  with  randomly  spaced  gaps  are  computed  and  their  effect 
on  observations  of  solar  oscillations  are  discussed.  It  is  found  that 
for  all  the  cases  considered,  observations  of  solar  oscillations  would 
not  be  severely  impacted  as  long  as  the  gap  structure  is  random  rather 

than  periodic. 

in  Solar  Seismology  from  Space  (eds.  R.K.  Ulrich  et  al.:  NASA/ 
Jet  Propulsion  Laboratory  Publ.  84-84),  pp.  255-262  (1984) 
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THE  EFFECTS  OF  IMAGE  MOTION  ON  THE  1- v  DIAGRAM 


Frank  Rill  # 
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Joint  Institute  for  Laboratory  Astrophysics, 
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ABSTRACT:  A  simple  two-dimensional  (x,t)  model  of  the  solar  oscillatory 
velocity  field  is  subjected  to  a  form  of  differential  image  motion. 
This  image  motion  is  meant  to  approximately  model  the  effect  of  the 
Earth's  atmosphere  on  observations  of  high  degree  solar  oscillations. 
The  distorted  velocity  field  is  analyzed  to  provide  the  apparent  fre¬ 
quencies  of  the  modes.  Comparison  of  the  results  with  the  frequencies 
obtained  from  the  undlstorted  case  shows  that  the  image  motion  can 
produce  a  discrepancy  of  as  much  as  12  pNz . 

^Operated  by  the  Association  of  Universities  for  Research  in  Astronomy, 
Inc.  under  contract  AST  78-17292  with  the  National  Science  Foundation 
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(ed.  D.O.  Gough;  ReideJ,  Dordrecht),  pp.  85-92  (1986) 
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INFLUENCE  OF  SPATIAL  FILTERING  ON 
POSSIBLE  ANISOTROPIES  IN  SOLAR  OSCILLATIONS 


Frank  Hill 
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Tucson,  Arizona  85726  USA 

Deborah  A.  Haber  and  Juri  Toomre 
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and 

Laurence  J.  November 

Air  Force  Geophysical  Laboratory 
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ABSTRACT.  We  have  used  full  disk  Doppler  observations  of  solar  oscillations 
to  compare  the  amplitudes  of  sectoral  modes  propagating  along  the  equator  with 
those  of  similar  modes  propagating  along  a  great  circle  aligned  with  the  poles.  We 
find  that  the  amplitudes  are  generally  not  equal  for  the  two  classes  of  modes,  but 
the  results  are  sensitive  to  analysis  procedures  attempting  to  isolate  the  different 
modes  of  oscillation.  Spatial  filtering  of  the  data  using  spherical  harmonics  sug¬ 
gests  that  greater  amplitudes  are  associated  with  “polar"  sectoral  modes  than  with 
“equatorial"  sectoral  modes. 


’n  Proc.  1AU  Symp.  123.  Advances  in  Helio-  and  Asteroseismoloqy 

\eds-  J.  Christensen-Dalsgaard  and  S.  Frandsen;  Reidel,  Dordrecht), 
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SOLAR  EQUATORIAL  ROTATION  RATE  INFERRED  FROM  INVERSION 
OF  FREQUENCY  SPLITTING  OF  HIGH-DEGREE  MODES 


Frank  Hill 
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Douglas  Cough 
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Theoretical  Physics,  University  of  Cambridge,  Cambridge  CB3  OHA,  England 

Jurl  Toomre  and  Deborah  A.  Haber 

Joint  Institute  for  Laboratory  Astrophysics,  and  Department  of  Astrophysica), 
Planetary  and  Atmospheric  Sciences,  University  of  Colorado, 

Boulder,  CO  80309-0440,  USA 

ABSTRACT  The  equatorial  rotation  rate  has  been  inferred  as  a  function  of  depth 
through  the  outer  16  Mm  of  the  Sun  from  observations  of  bigh-degree  five-minute  oscilla¬ 
tions  An  optimal  averaging  inversion  procedure  due  to  Backus  L  Gilbert  (1970)  has  been 
applied  to  frequency  splittings  measured  from  power  spectra  obtained  using  Doppler  data 
spanning  three  and  five  consecutive  days.  The  resulting  rotation  curves  have  proven  to 
be  much  more  stable  than  the  curves  obtained  from  data  sets  of  single  days  The  results 
imply  that  the  solar  rotation  rate  increases  with  depth  by  0  023  pHi  reaching  a  maximum 
at  about  2  Mm  below  the  surface,  then  decreases  by  0.037  pHx  down  to  16  Mm 


v’o'.yvv.v.'Avv.v 


Astrophys.  Jt,  to  be  submitted  (1987) 

PAPER  M. 

EFFECTS  OF  MODE  BEATING  ON  DETERMINATION 
OF  SOLAR  OSCILLATION  FREQUENCIES 
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and 
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ABSTRACT.  Numerical  simulations  have  been  used  to  assess  effects  of  mode 
beating,  or  interference  between  unresolved  modes  of  solar  oscillation,  upon  the  de¬ 
termination  of  their  frequencies.  Synthetic  data  sets  of  oscillations  are  constructed 
for  modes  with  spherical  harmonic  degrees  t  ranging  from  50  to  400.  The  isola¬ 
tion  of  sectoral  modes  is  sought  by  using  a  narrow  rectangular  observing  window 
aligned  with  the  equator  and  averaging  the  data  in  longitude  before  carrying  out 
the  Fourier  transforms  to  determine  power  spectra.  The  simulations  suggest  mode 
beating  can  produce  differences  of  about  15  pHz  between  sectoral  mode  frequencies 
measured  from  ridge  centroids  in  the  power  spectra  and  the  actual  input  frequen¬ 
cies.  Variations  in  mode  amplitudes  and  phases  from  one  realization  to  another  in 
these  simulations  can  lead  to  differences  of  about  10  pHz  when  the  measured  fre¬ 
quencies  are  compared.  These  errors  are  a  significant  source  of  noise  for  inversion 
of  helioseismic  data. 


(  V,*  *,l  .  ‘  c  v.t'v.s,  ’  *  .4' 1  ■* 


TW  ASTKOmmc*!.  Kx-»sw_a3  S79-I86.  !9M  Aufuit  15 
f  mi  TIL.  ITITT . . PlKWllClA 


Astrophys.  J.»  283,  879-886  (1984) 


PAPER  N, 


FREQUENT  ULTRAVIOLET  BR1GHTEN1NGS  OBSERVED  IN  A  SOLAR  ACTIVE  REGION 

WITH  SOLAR  MAXIMUM  MISSION 

Jasov  G  Porte*.  Ju*i  Toomke.  and  Katharine  B  Gebbie 
*».  fc,  L»b©f»iot)  w»».' 
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ABSTRACT 

rkWrv.tmr.*  in  the  ultraviolet  of  sites  of  enhanced  intensity  within  an  active  region  on  the  Sun  reveal  fre- 
OcTs  wmphng  in" ime'S^he UUwioVet' instrument  on  the  Solar 
Maximum  telectitg^the'brightes^sUe'aratj^ventime 

It  appears  that  the  modulations  in  intensity  are  found  to  correlate  with  the  jitter  tn  the  satel- 

Jeauig  due  to  magnetic  field  reconnection  within  an  active  region  is  proceeding  almost  stochastically .  events 
involving  only  a  modest  release  of  energy  occur  the  most  frequently. 

Subject  headings:  Sun:  activity  —  Sun:  spectra  -  ultraviolet:  general  —  ultraviolet,  spectra 
'  ||LA  it  operated  jointl)  b)  the  Nitionel  Bureivi  of  Standards  and  the 
L’mvenit)  of  Colorado 

in  Proc.  IAU  Symp.  123,  Advances  in  Helio-  and  Asterosei sinology 


(eds.  J.  Christensen-Dalsgaard  and  S.  Frandsen;  Reidel,  Dordrecht), 
in  press  (1987) 
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RESPONSE  OF  THE  SOLAR  FIVE-MINUTE  OSCILLATIONS 
TO  A  MAJOR  FLARE 


Deborah  A.  Haber  and  Juri  Toomre 

Joint  Institute  for  Laboratory  Astrophysics,  and  Department  of  Astrophysical, 
Planetary  and  Atmospheric  Sciences,  University  of  Colorado, 

Boulder,  CO  80309-0440,  USA 

Frank  Hill 

National  Solar  Observatory,  National  Optical  Astronomy  Observatories, 
Tiicaon,  AZ  85726-6732,  USA 

ABSTRACT  —  Solar  five-minute  oscillations  of  intermediate-degree  l  were  observed  both 
before  and  after  a  very  strong  white-light  flare.  Intensity  images  of  the  full  Sun  taken  on 
two  sides  of  the  Fe  I  A  5576  spectral  line  were  recorded  on  film,  digitized  with  8*  spatial 
resolution,  and  then  converted  into  Doppler  velocities.  The  data  were  projected  onto  both 
equatorial  and  polar  sectoral  modes  and  Fourier  transformed  in  time.  Comparing  the 
resulting  power  spectra,  we  find  a  substantial  increase  in  power  in  the  ps  ridge  of  the 
equatorial  modes  on  the  day  after  the  flare;  such  an  increase  may  be  a  consequence  of  the 
solar  flare.  When  data  from  all  the  ridges  are  considered,  there  is  an  average  increase  in 
power  of  only  a  few  percent  the  day  after  the  flare.  This  overall  increase  is  probably  not 
significant  due  to  uncertainties  from  effects  of  the  beating  of  unresolved  modes. 
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Penetrative  cellular  convection  in  a  stratified  atmosphere 
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Sammary.  Penetrative  convection  has  been  investigated  within  a 
simple  compressible  model  consisting  of  three  layers  of  diflering 
stratification  prior  to  the  onset  of  convection.  The  middle  one  is 
a  cooveaively  unstable  polytrope  bounded  above  and  below  by 
two  stably  stratified  polytropes.  Nonlinear  anelastic  modal 
equations  are  used  to  describe  these  convective  motions.  Penetra¬ 
tive  boundary  conditions  are  imposed  at  the  top  and  bottom  of 
the  computational  domain.  One-  and  two-mode  steady  solutions 
with  bexagonsl  planforms  have  been  studied  for  Rayleigh  num¬ 
bers  up  to  about  10*  timet  critical,  and  for  a  range  of  Prandtl 
numbers,  horizontal  wavenumbers  and  stratifications.  These 
show  that:  1)  Penetration  into  the  lower  stable  layer  by  down¬ 
ward-directed  plumes  is  substantially  larger  in  a  stratified 
medium  than  in  a  Boussinesq  fluid,  and  produces  an  extended 
region  of  adiabatic  stratification.  2)  Overshooting  into  the  upper 
stable  layer  can  be  sharply  diminished  by  a  reversal  of  the 
buoyancy  occurring  in  tbe  upper  part  of  the  unstable  zone  if  the 
stratification  is  strong  enough.  These  efiects  are  all  related  to 
large  pressure  perturbations  that  arise  in  such  convective  flows. 
The  strong  asymmetry  between  upward  and  downward  penetra¬ 
tion  in  compressible  media  bas  major  implications  for  the  mixing 
of  stable  regions  above  and  below  stellar  convection  zones. 

Bey  words:  stellar  structure  -  convection  zones  -  astrophysical 
fluid  dynamics 
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STRONG  DOWNWARD  PLUMES  RESULTING  FROM  COMPRESSIBILITY  IN 
NONLINEAR  CONVECTION  AND  THEIR  COUPLING  TO  GRAVITY  WAVES 


JURI  TOOMRE,  NEAL  E.  HURLBURT,  and  JOSEP  M.  MASSAGUER  1 

Joint  Institute  for  Laboratory  Astrophysics,  and 
Department  of  Astrophysical ,  Planetary  and  Atmospheric 
Sciences,  University  of  Colorado,  Boulder  CO  80309 


ABSTRACT 

Two-dimensional  numerical  simulations  are  used  to  model  fully 
compressible  nonlinear  convection  spanning  multiple  scale  heights  with¬ 
in  a  stellar  envelope.  The  resulting  cellular  flows  display  prominent 
downward -directed  plumes  surrounded  by  broader  regions  of  upflow.  Such 
asymmetry  arises  because  pressure  fluctuations  accentuate  buoyancy 
driving  In  the  concentrated  plumes  and  can  lead  to  buoyancy  braking  in 
the  surrounding  ascending  flows.  When  such  convection  is  allowed  to 
overshoot  into  regions  of  stable  stratification  above  and  below  the 
unstable  zone,  then  the  dynamics  of  penetration  is  dominated  by  the 
downward  plumes.  This  leads  to  the  excitation  of  a  broad  spectrum  of 
internal  gravity  waves  in  the  lower  stable  zone,  and  these  waves  in 
turn  can  produce  a  significant  modulation  in  the  amplitude  of  the 
convection  through  their  feedback  upon  the  plumes. 


*Also  at  E.T.S.  Ingenieros  de  Caminos,  Universidad  Polytecnica  de 
Barcelona,  Barcelona,  Spain. 
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ON  THE  DETERMINATION  OF  THE  LIFETIME  OF  VERTICAL  VELOCITY 
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ABSTRACT 

Observational  studies  of  the  vertical  velocities  of  mesogranula- 
tion  and  supergranulation  provide  conflicting  results  for  the  lifetimes 
of  these  patterns  when  analyzed  by  two  different  methods.  Visual  in¬ 
spection  of  the  velocity  images  suggests  that  mesogranulation  has  a 
lifetime  in  excess  of  2  hours,  while  cross-correlation  methods  imply  a 
lifetime  of  only  about  40  min.  For  supergranulation,  the  correlation 
technique  yields  a  lifetime  of  2.7  hours,  far  short  of  the  24  hours 
found  by  many  other  studies  considering  network  structures  or  horizon¬ 
tal  velocities.  These  discrepancies  may  be  due  to  temporal  changes  in 
the  apparent  velocity  zero,  by  differential  distortions  in  the  images 
due  to  large-scale  flows,  or  by  evolution  in  the  flows  controlled  by 
the  magnetic  fields. 


1  Operated  by  the  Association  of  Universities  for  Research  in 
Astronomy,  Inc.  under  contract  AST  78-17292  with  the  National  Science 
Foundation 

2  NAS/NRC  Resident  Research  Associate 
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TWO-DIMENSIONAL  COMPRESSIBLE  CONVECTION  EXTENDING  OVER 
MULTIPLE  SCALE  HEIGHTS 
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ABSTRACT 

Two-dimensional  simulations  are  used  to  study  fully  compressible  thermal  convection  spanning  multiple 
density  scale  heights  typical  of  a  stellar  envelope.  The  fluid  is  assumed  to  be  a  perfect  gas  with  constant 
thermal  conductivity  and  dynamic  viscosity.  The  unstably  stratified  layer  is  thus  represented  by  a  polytropic 
index  typically  taken  to  be  n  *  1  (with  a  *  3/2  equivalent  to  an  adiabatic  stratification).  The  mean  density 
ratio  (bottom  to  top  of  the  layer)  ranges  from  about  I  to  21,  with  Rayleigh  numbers  up  to  about  1000  times 
critical  These  highly  nonlinear  flows  are  studied  with  a  two-dimensional  numerical  scheme  based  on  a  modi¬ 
fied  two-step  Lax-WendrofT  method  The  convective  motions  in  these  simulations  span  the  full  height  of  the 
unstable  layer,  with  no  tendency  to  form  a  succession  of  rolls  in  the  vertical  as  has  been  assumed  in  mixing- 
length  treatments  of  convection  Further,  the  flows  remain  subsonic  because  the  center  of  the  cell  shifts  toward 
the  bottom  of  the  layer  as  the  density  stratification  is  strengthened.  The  flows  then  display  prominent 
downward-directed  plumes  surrounded  by  broader  regions  of  upflow  The  flow  asymmetry  leads  to  a  kinetic 
energy  flux  which  is  directed  downward,  whereas  the  enthalpy  or  convective  flux  is  upward  Such  asymmetry 
comes  about  because  pressure  fluctuations  accentuate  buoyancy  driving  in  the  downward  plumes  and  can  lead 
to  buoyancy  braking  in  the  surrounding  ascending  flows.  Compressions!  work  is  significant  in  the  overall 
energetics. 

Subject  headings  convection  —  hydrodynamics  —  stars :  atmospheres 
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ABSTRACT 

Penetrative  convection  spanning  multiple  scale  heights  is  studied  within  a  simple  stellar  envelope  consisting 
of  three  layers  a  convectively  unstable  middle  layer  bounded  above  and  below  by  stably  stratified  polytropes 
Two-dimensional  numerical  simulations  are  used  to  investigate  the  fully  compressible  nonlinear  motions  that 
ensue  The  cellular  flows  display  prominent  downward-directed  plumes  surrounded  by  broader  regions  of 
upflow  Such  asymmetry  arises  because  pressure  fluctuations  accentuate  buoyancy  driving  in  the  concentrated 
plumes  and  can  even  lead  to  weak  buoyancy  braking  in  the  surrounding  ascending  flows  As  the  plumes 
plunge  downward  into  a  region  of  stable  stratification,  they  serve  to  excite  a  broad  spectrum  of  interna! 
gravity  waves  there  The  induced  waves  are  not  passive,  for  they  feed  back  upon  the  plumes  by  deflecting 
them  sideways,  thereby  modulating  the  amplitude  of  the  convection  in  time  even  in  the  unstable  layer  The 
penetrative  motions  that  billow  upward  into  the  upper  stable  zone  are  distinctly  weaker,  and  they  cascade 
back  downward  toward  the  unstable  zone  over  a  broad  horizontal  scale  The  strong  excitation  of  grant; 
waves  by  the  convection  has  implications  for  gradual  mixing  deep  within  a  star 
Subject  headings:  convection  —  hydrodynamics  —  stars  interiors 
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ABSTRACT 

Two-dimensional  numerical  simulations  are  used  to  study  fully  compressible  con¬ 
vection  in  the  presence  of  an  imposed  magnetic  field.  Highly  nonlinear  flows  are  con¬ 
sidered  that  span  multiple  density  scale  heights.  The  convection  tends  to  sweep  the 
initially  uniform  vertical  magnetic  field  into  concentrated  flux  sheets  with  significant 
magnetic  pressures.  These  flux  sheets  are  partially  evacuated,  and  effects  of  buoyancy 
and  Lorentz  forces  there  can  serve  to  suppress  motions.  The  flux  sheets  can  be  sur¬ 
rounded  by  a  sheath  of  descending  flov  If  the  imposed  magnetic  field  is  sufficiently 
strong,  the  convection  displays  a  variety  of  oscillations.  The  unstably  stratified  fluid 
layer  has  an  initial  density  ratio  (bottom  to  top  of  layer)  of  11.  Surveys  of  solutions  at 
fixed  Rayleigh  number  sample  Chandrasekhar  numbers  from  1  to  1000  and  magnetic 
Prandtl  numbers  from  1/16  to  1.  These  nonlinear  simulations  utilize  a  two-dimensional 
numerical  scheme  based  on  a  modified  two-step  Lax-Wendroff  method. 
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Abstract 

Two-dimensional  numerical  simulations  are  used  to  investigate  how  fully  com¬ 
pressible  nonlinear  convection  penetrates  into  and  mixes  a  stably  stratified  zone 
beneath  a  stellar  convection  sone.  The  dynamics  of  penetration  is  dominated  by 
downward  plumes  which  can  extend  far  into  the  stable  material  and  which  lead 
to  the  excitation  of  a  broad  spectrum  of  internal  gravity  waves  in  the  lower  stable 
sone.  The  depth  of  penetration  is  controlled  by  a  balance  between  the  kinetic  en¬ 
ergy  carried  into  the  stable  layer  by  the  plumes  to  and  the  buoyancy  braking  they 
experience  there.  The  motions  mix  the  chemical  stratification  in  the  penetrative 
layer  in  a  few  convective  overturning  times.  The  model  fluid  is  taken  to  be  a  per¬ 
fect  gas  with  a  constant  dynamic  viscosity.  The  computational  domain  is  divided 
into  regions  of  stable  and  unstable  thermal  stratifications  by  varying  the  thermal 
conductivity  with  depth.  The  mean  density  ratio  (bottom  to  top  of  the  layer)  is 
initially  6  across  the  unstable  sone  and  114  across  the  entire  computational  domain, 
with  Rayleigh  numbers  in  the  unstable  layer  taken  to  be  about  100  times  critical. 
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Recent  calculations  suggest  that  the  bulk  of  the  solar  magnetic  Held  may  be  stored  in 
a  thin  convectively  stable  region  situated  between  the  convection  2one  proper  and  the 
radiative  zone.  Determining  the  stability  properties  of  such  a  field  is  therefore  impor¬ 
tant  with  implications  for  both  the  generation  and  escape  of  magnetic  flux  The 
magneto-Boussinesq  equations  are  used  to  perform  a  linear  stability  anal) sis  of  a 
static  plane  layer.  Several  new  instability  mechanisms  are  revealed  showing  previous 
ideas  concerning  magnetic  buoyancy  instabilities  to  be  over  simplified.  The  most 
important  result  is  that  instability  may  occur  even  for  fields  which  increase  with 
height.  Detailed  results  are  presented  for  2  and  3  dimensional  motions  for  a  weakly 
stratified  magnetic  field  together  with  a  simple  calculation  for  the  2  dimensional 
instability  of  a  strongly  varying  field. 
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Recent  calculations  suggest  that  the  bulk  of  the  solar  toroidal  field  may  be  stored  in  a 
thin,  convectively  stable  region  situated  between  the  convection  zone  proper  and  the 
radiative  zone.  Determining  the  stability  properties  of  such  a  field  is  therefore 
important  with  implications  for  both  the  generation  and  escape  of  magnetic  flux.  The 
plane  layer,  linear  stability  analysis  of  Hughes  (1985)  is  extended  to  incorporate  the 
effects  of  uniform  rotation.  Detailed  studies  are  made  of  interchange,  or 
“axisymmetric"  modes  and  of  undular,  or  wavelike,  motions,  considering  modes  of 
both  low  and  high  frequency.  The  force  due  to  rotation  acts  to  constrain  the  fluid 
motions,  a  feature  which  is  strongly  stabilizing  for  direct  modes,  but  can,  in  certain 
circumstances,  be  destabilizing  for  oscillatory  motions. 

For  the  interchange  modes  we  show  that  the  instability  discussed  at  length  by 
Hughes  (1985),  driven  by  fields  increasing  with  height,  is  still  present  and  indeed  may 
be  enhanced  by  rotational  effects.  We  also  study  the  more  conventional  instabilities, 
discussing  the  transformation  between  direct  and  oscillatory  modes  and  considering  in 
detail  some  peculiar  properties  of  the  oscillatory  instabilities. 

The  more  relevant  instabilities  in  an  astrophysical  context  are  likely  to  be  undular 
modes.  Previous  studies  of  low  frequency  modes  driven  by  top  heavy  field  gradients 
are  extended  to  consider  modes  of  various  frequencies  for  a  wide  range  of  parameter 
values.  Of  particular  interest  is  the  occurrence  of  two  distinct  modes  of  instability  for 
bottom  heavy  field  gradients.  We  also  exhibit  some  of  the  peculiar  stability 
boundaries  which  can  result  when  none  of  the  competing  influences  in  the  problem  is 
dominant. 
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Magnetic  buoyancy  instabilities  may  be  a  significant  cause  of  flu*  loss  from  the  solar 
interior  This  paper  extends  the  hnear  theory  for  such  instabilities,  which  has  been 
widely  studied,  into  the  weakly  non-linear  regime  to  obtain  solutions  valid  for  small 
but  finite-amplitude  motions  We  consider  both  static  and  routing  equilibrium  sutes 
but.  for  simplicity,  focus  attention  only  on  the  two-dimensional  interchange  modes 
The  method  used  in  obuinmg  the  solutions  is  due  to  Coullei  and  Spiegel  (1983)— a 
summary  of  the  method,  together  with  its  advantages,  is  contained  in  Section  2. 

For  the  static  equilibrium  the  bifurcation  structure  is  similar  to  that  of  thenno- 
haline  convection  and  magnetoconvection  and  »  fairly  straightforward  to  analyte  The 
inclusion  of  rotation  yields  more  complicated  non-linear  behaviour  with  the  ap¬ 
pearance  of  more  varied  unfoldings  as  well  as  tome  new  bifurcations  One  of  the  key 
results  to  emerge  is  that  for  both  static  and  rotating  basic  turn  the  oscillatory  modes 
destabilised  by  'bottom-heavy*  field  gradients  (magnetic  field  increasing  with  height) 
art  often  unstable  to  finite-amplitude  disturbances  although  stable  to  infinitesimal 
perturbations  A  detailed  discussion  of  all  the  finite-amplitude  behaviour  is  given  in 
Sections  3  and  4 

KEY  WORDS  Magnetic  buoyancy  instability,  bifurcation,  weakly  non-linear  regime, 
equilibrium  states 
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ABSTRACT 

Although  the  undular  instabilities  of  a  stratified  horizontal  magnetic  field  have 
been  studied  in  a  number  of  contexts  we  believe  that  the  physical  mechanism  responsible 
for  the  instability  has  not  been  satisfactorily  explained.  In  this  paper  we  present  a  new 
explanation  of  why  these  instabilities  occur,  considering  in  detail  the  quite  different  cases 
of  two-dimensional  and  three-dimensional  motions. 
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Abstract 

The  connection  between  the  effective  convection  velocities  for  scalar  and  vector 
fields  due  to  the  action  of  turbulence  is  discussed.  An  explicit  relation  between  the  two 
is  calculated  for  the  important  special  case  of  two-dimensional  Sows  and  fields  and  it  is 
shown  that  both  velocities  are  sero  for  homogeneous  turbulence.  The  calculation  leads  to 
new  insights  on  the  r61e  of  symmetry  in  determining  the  advection  velocity  for  the  case  of 
a  vector  field. 
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ABSTRACT 

In  order  that  magnetic  flux  be  confined  within  the  eolar  interior  for  time*  comparable  to 
the  tolar  cycle  period  it  ha*  been  suggested  that  the  bulk  of  the  tolar  toroidal  field  ia  stored 
in  the  convectivcly  *tabl*  overshoot  region  situated  beneath  the  convection  tone  proper.  Such  a 
magnetic  field  though  is  still  buoyant  and  is  therefore  subject  to  Rayleigh-Taylor  type  instabilities, 
la  this  paper  we  consider  the  model  problem  of  an  isolated  region  of  magnetic  field  embedded  in  a 
convectivcly  stable  atmosphere.  The  fully  nonlinear  evolution  of  the  two-dimensional  interchange 
modes  is  studied,  thereby  shedding  some  light  on  one  of  the  processes  responsible  for  the  escape 
of  flux  from  the  solar  interior. 
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The  partial  differential  equations  governing  two  dimensional  thericosolutal  convec¬ 
tion  in  a  Boutsinesq  fluid  with  free  boundary  conditions  have  been  solved  numerically 
m  a  regime  where  oscillatory  solutions  can  be  found  A  systematic  study  of  the 
transition  from  nonlinear  periodic  oscillations  Vo  temporal  chaos  has  revealed 
sequences  of  period-doubling  bifurcations  Overstsbility  occurs  if  the  ratio  of  the 
aolutal  to  the  thermal  diffusivity  r<  I  and  the  soluul  Rayleigh  number  Rs  is 
sufficiently  large  Solutions  have  been  obtained  for  two  representative  values  of  r 
For  r  -  0  316,  Rs  «  10*.  symmetrical  oscillations  undergo  a  bifurcation  to  asym¬ 
metry.  followed  by  a  cascade  of  period-doubling  bifurcations  leading  to  aprnodiritv. 
as  the  thermal  Rayleigh  number  RT  is  increased  At  higher  values  of  RT.  the 
bifurcation  sequence  is  repeated  in  reverse,  restoring  simple  periodic  solutions  As  RT 
it  further  increased  more  period-doubling  cascades  followed  by  chaos,  can  be 
identified  Within  the  chaotic  regions  there  are  narrow  periodic  windows  and 
multiple  branches  of  oscillatory  solutions  coexist  Eventually  the  oscillatory  branch 
ends  and  only  steady  solutions  can  be  found  The  development  of  chaos  has  been 
investigated  for  r  -  0  I  by  varying  Rr  for  several  different  values  of  Rs  When  Rs 
it  sufficiently  small  then-  an  periodic  solutions  w  hose  period  becomes  infinite  at  the 
end  of  the  oscillatory  branch  As  R s  is  incnaaed  chaos  appears  in  the  neighbourhood 
of  these  hfteroclinic  orbits  At  higher  v  slues  of  Rs  chaos  is  found  for  a  broader  range 
in  Rt  A  truncated  fifth  order  model  suggests  thst  the  spfversncr  of  chaos  is 
aasocisted  with  hfteroclinic  bifurcations 


